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ABSTRACT 

We examine the central-galaxy luminosity - host-halo mass relation for 54 Brightest 
Group Galaxies (BGGs) and 92 Brightest Cluster Galaxies (BCGs) at z < 0.1 and 
present the first measurement of this relation for a sample of known BCGs at 0.1 < 



z < 0.8 [z ~ 0.3). At z < 0.1 we find Lk (x M. 



0.24±0.08 
200 



for the BCGs and the early- 



type BGGs in groups with extended X-ray emission and Lk oc M. 



0.11±0.10 
200 



for the 



BCGs alone. At 0.1 < z < 0.8 we find Lk oc M°oo^^°-". We conclude that there is no 
evidence for evolution in this relationship between z < 0.1 and z < 0.8: BCG growth 
appears to still be limited by the timescale for dynamical friction at these earlier times, 
not proceeding according to the predictions of current semi-analytic models. 

Key words: Galaxies: clusters: general - galaxies: elliptical and lenticular, cD - 
galaxies: evolution - galaxies: formation 
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1 INTRODUCTION 

Brightest cluster galaxies (BCGs) are the most mas- 
sive galaxies known in the Universe and are observed 
to have remarkably small dispersion in their absolute 
mag nitudes compared to other early-type cluster galaxie s 
(e.g. lAragon-Salamanca et ai]|l998l : ICoUins fc Mannill998h . 



Their position at the centre of clusters suggests their evolu- 
tion and their uni que propert i es are linked to their specia l 
environment (e.g. lEdgd Il99il : Ivon der Linden et al]|2007l ). 
However, it is still not clear exactly which processes drive 
their evolution. 

Of central importance in this context is the mass assem- 
bly history of these most massive galaxies and how that de- 
pends on their environment. Semi-analytic models based on 
merger trees from the Millennium N-body simulation predict 
that, while the stars in BCGs form at high-redshift, the mass 
of the g alaxy only a ssembles recently (doubling in mass since 
z ~ 1; Ide LuciaTfc Blaizot 2007 ), with the most massive 
galaxies in the most massive clusters having undergone the 
most mass assembly. Observations of BCGs agree that their 
prop erties are consistent with having undergone me rgers 
(e.g. iBovlan-Kolchin et al.l [20061 : iBernardi etall |2007| ) and 
BCGs in the most massive clusters do appear to have un- 
dergone more stellar mass evolution than those in less mas- 
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sive clusters, as shown by their surface brightness profiles 
(Brough et al. 2005). However, that evolution must have 
predominantly occurred before z ^ 1 as shown by the un i- 
formity in their absolute mag nitudes llBrough et al.l 120021 ') , 
and their steep metallicity gradients iBrough et al.ll2007l ). 

An important clue in this context lies in the central 
galaxy luminosity-host-halo mass (LcMh) relation, which 
shows a scaling b e tween BCG mass and that of its host 
halo. iLin fc Mohrl ||2004') found that Lk oc 

oXu ■ 



0.26±0.04 



for 



93 BCGs at z < 0.1, iPopesso et al.1 (|2007l ) found Lr oc 
^^200^"^° °" for 217 BCGs at 0.0 1 < z < 0.25 in the Sloan 
Digital Sky Survey (SDSS), and [Hansen et al.l l|2007l ) found 
{Li) oc M°og°±'' °i for the SDSS MaxBCG cluster sample 
(0.1 < z < 0.3). Studies of the LcMh-relation have not ven- 
tured beyond z ~ 0.3, so it is not clear whether this rela- 
tionship evolves with cosmic time. 

The LcMh-relation is also a natural output from anal- 
yses of cosmological simulatio ns: Halo Occupation Dis- 
tribution function s (HOD; e.g. Ivan den Bosch et al. I l2005l : 
IZheng et al.l |2007| ) and Conditional Luminosity Functions 
(CLF; e.g. Cooray & Milosavljevic 2005; henceforth CM05). 
These statistically assign numbers of galaxies or the lumi- 
nosity distribution of galaxies in a given dark matter halo, 
distinguishing between central and satellite galaxies and en- 
abling a direct comparison with observations. A related ap- 
proach is the application of a direct relationship between 
the mass of a dark matter halo or sub-halo and the ob- 
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served galaxy luminosity (e.g. IVale Ostrikerl |2006| . l2007l: 
henceforth VO06, VO07). CM05 predict (L) oc M^Jl,-^ above 
halo masses of 4 x 10^^h~^ Mq as the timescale on which 
dynamical friction causes satellite galaxies to fall in to the 
centre exceeds the age of the host system halo. VO06 pre- 
dict {Lk) oc M?of as a result of the hierarchical formation 
of structure in the Universe, with a prescription for sub-halo 
mass loss. We interpret the LcMh-relation in this paradigm: 
that it results from the BCG growing in step with its host 

halO; 

lYang et all (120051 ') analysed the HOD of galaxy groups 
in the 2-degree Field Galaxy Redshift Survey ( 0.01 < z < 
0.2) a nd found {Lcm,bj) oc M^'^^ . More recentlv. lYang et al.l 
|2003) determined the CLF of groups (0.01 < z < 0.2) in 
the SPSS and foun d {Lcen,r) oc Af^'^^. At higher redshifts, 
IZheng et all (|2007l ) analysed the HOD using the two-point 
correlation function of the DEEP2 galaxy redshift survey 
[z ~ 1) in comparison to the SDSS {z ~ 0) and determined 
the LcMh-relation and its evolution. They found the VO06 
model to be a good fit to their data at z ~ a nd z ^ 1 sug- 

Jesting that this relationship has not evolved. IWhite et al.l 
200% found (Lcen.s) oc M^-^^ at 2 ~ 0.5 from their HOD 
analysis of ~ 2L* red galaxies in the NOAO Deep Wide 
Field Survey. These studies, however, do not directly iden- 
tify central galaxies. 

In this Letter, we examine the central-galaxy luminos- 
ity - host halo mass relation for a sample of known BCGs 
in groups and clusters at 2: < 0.1 and, for the first time, 
the relationship for known BCGs in clusters at redshifts 
0.1 < z < 0.8. This will enable us to directly determine 
the evolution in this relationship for known BCGs and to 
compare to the HOD ana lyses at these redshi f ts. It will also 
enable us to examine the lde Lucia fc BlaizotI (|2007l ) predic- 
tion of these galaxies having doubled in mass since z ~ 1. 

In Section [2] we introduce our sample and then present 
our comparison of BGG and BCG luminosities with host 
system mass in Section We discuss the implications of our 
results in Section |4] and draw our conclusions in Section (5] 
Throughout this paper we assume Ho = 70 km s"'^ Mpc~^, 
Qm = 0.3 and = 0.7. 



2 DATA 

The BCG sample is from lBrough et all (|2002| ) and contains 
92 BCGs at z < 0.1 and 63 at 0.1 < z < 0.8 {z ~ 0.3). 
These are selected as the brightest galaxies closest to the 
X-ray centroid of their host cluster and are all early-type 

galaxies^ 

InHrougheiaD (|2002l ). if -band magnitudes were mea- 
sured in an aperture of radius 12.5/i~^ kpc. The apparent 
magnitudes were corrected to absolute magnitudes using K- 
and passive evolution corrections (Yoshii & Takahara 1988) 
produced by the GISSEL9 6 stellar population synthesis code 
IIBruzual fc Charioting ). The galaxies were assumed to be 
10 Gyrs old, to have formed in an instantaneous burst, and to 
have evolved passively since 2 ~ 2. The assumptions of age 
and redshift of formation have negligible effect on the correc- 
tion value. The magnitudes were also corrected for galactic 
ab sorption (typical valu es Ak ^ 0.02 mag) using the maps 
of lSchlegel et al. I l|l998l ). The uncertainties in these magni- 
tudes are 0.04 mag. 



Cluster X-ray luminosity is directly proportional to the 
square of the electron density of the intra-cluster medium 
and provides a measure of the host system mass. X-ray 
luminosities for the low-redshift clusters were measured 
from the ROSAT All Sky Survey (RAS S), paired wit h co- 
ordinates from AbeU et al. (1989) and iLvnanJ Il999l . The 
X-ray luminosities for the high-redshift clusters were ob- 
tained from the E instein Medium Sensitivity Survey (EMSS; 
iGioia fc Luppinol [l994l ) and Serend ipitous High-Red shift 
Archival Cluster j^Q g^T catalogues l|Burke et al.ll2003l ). In 
iBrough et al] (|2002h the RASS passband (0.1-2.4 kev) data 
were transformed to the EMSS passband (0.3 — 3.5 keV). 
Here we use the same transformation to place all these values 
in the RASS passband, i.e. Lx (0.1 -2.4 keV)= Ljf(0.3-3.5 
keV)/1.08. The host clusters cover a wide range in X-ray lu- 
minosity: 5 X 10**^ < Lx < 5 X 10'*'* {h^^ erg s~^) at redshifts 
z < 0.1 and 6 x 10*^ < Lx < 1 x 10''^ (ft^^ erg s~^) at red- 
shifts 0.1 < z < 0.8. 

The Group Evolution Multiwavelength S urvey (GEMS; 
lOsmond fc Ponmanlliooi : iForbes et aLll2006l ) is a heteroge- 
neous sample of 60 groups selected from optically-selected 
group catalogues observed by the ROSAT X-ray satellite to 
enable X-ray classification. Removing the 6 groups for which 
lOsmond fc PonmanI l|2004l ) were unable to find at least 4 
galaxies associated with the original optical position, leaves 
a sample of 54 groups with 2 x 10*° < Lx < 4 x 10*^ (/i"^ 
erg s~^). Of these 54 groups, 35 have extended intra-group 
X-ray emission, 13 have X-ray emission only associated with 
the central galaxy and 6 are undetected (at 3 a above back- 
ground) in X-rays. lOsmond fc PonmanI l|2004h also identified 
the BGG and published the galaxies' morphologies. 

All the GEMS groups are covered by the 2 Mi- 
cron Ah Sky Surv ey Extended Source Catalogue (2MASS; 
I jarrett et ahlbOOOl ). Our photometric reductions w ere per- 
forme d to be as consistent as possible with Broug h et al.l 
(j2002l ). Ah 54 groups are at redshifts z < 0.03 {z ~ 0.01) 
and for a large proportion of the sample the circular aperture 
photometry in 2MASS does not extend to large enough radii 
to measure an aperture of radius 12.5h~^ kpc (the maximum 
possible radius 70" — 14h~^ kpc at 2 ~ 0.01). The aperture 
magnitudes for this sample were therefore measured manu- 
ally from the 2 MASS scans and corr ected to absolute magni- 
tudes following iBrough et aLll|2002l ). They have photometric 
uncertainties of ~ 0.04 mag. 

The group X-ray luminosities are taken from the GEMS 
group catalogue. These are bolometric luminosities, extrap- 
olated to the radius corresponding to an overdensity of 
500 times the critical density (rsp p), with uncertainties of 
logLx ~ 0.05. iBrough et al] ^20071 ) showed that the GEMS 
group X-ray luminosities are consistent with those measured 
by the ROSAT- ESO Flux Limited X-ra y galaxy cluster sur- 
vey (REFLEX; iBohringer et all |2004| ) extrapolated to 12 
times the core radius of the cluster, with an uncertainty 
of logLx = 0.05. We repeated our low-redshift group and 
cluster analysis using REFLEX X-ray luminosities for the 
40 clusters that are in common with this work and found 
that our results are unchanged. We therefore present this 
work using the full low-redshift sample of 92 clusters. 

We convert the group X-ray luminosities to the mass 
within an overdensity of 200 times the critical density, M200 , 
using the Lx (rsoo , Bol) — Mi^nrf r e lation measured for 
15 GEMS groups from IBrough et al.l (|2006l ). We then use 
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the Lx(0.1 - 2.4 keV)-M\-^^ ^°-°^ relation for clusters from 
iReiprich fc Bohringei] l|200a i for our clusters. We estimated 
the uncertainties in our mass measurements by perform- 
ing a Monte-Carlo analysis of how the uncertainties in the 
Lx measurements propagate through to our mass estimates 
given the uncertainty in the slope of the — M200 relation 
used for our conversion. These were found to be of order 
0.08 in log M200 fo r the groups and 0.02 for the clusters. 
lEttori et all \20Q4 \ show that the iReiprich fc Bohringei 
l|2002h relationship does not evolve significantly over the 
redshift range analysed here, we therefore use the cluster 
Lx(0.1 — 2.4 keV)— M200 relationship for our higher redshift 
sample. 

2.1 Aperture Magnitudes 

In their BCG studv. lLin fc Mohij (j2004h make the case that 
aperture magnitudes are not sensitive to the mass growth of 
BCGs and instead use 2MASS Tf-band isophotal magnitudes 
(measured within the 20th magnitude isophote). We test 
whether this affects our data in several different ways: 

We compared t he apparent magnitu des for the 24 galax- 
ies in common with iLin fc Mohil l|2004l ) and found an offset 
of Arriif = 0.35 ±0.03 mag resulting from the use of the dif- 
ferent measurements (Figure [T|. This offset is not correlated 
with the magnitudes or with cluster mass. 

We re-measure the magnitudes in a larger aperture 
of 25/i~^kpc for 4 BCGs at the same mean luminosity in 
high-mass (M ~ 4 X lO"/i"^M0), and in low -mass clusters 
we were subject to aperture effects 
then the magnitudes of the galaxies in the low-mass clusters 
would increase less than those in the high-mass clusters. 
However, the increases are consistent: AMif(low mass) = 
0.39 ± 0.02, AMif (high mass) = 0.44 ± 0.04. 

Finally, we use the fl -band surface brightness profiles 
from iGraham et all l|l996l ) of 24 galaxies in common with 
this work. We then predicted (assuming R~ K = 2.6) aper- 
ture luminosities within r = 12.5h~^ kpc and isophotal lu- 
minosities within K — 20 magarcsec"^. The lower panel of 
Figure [1] shows that aperture luminosities follow the same 
relationship as isophotal luminosities with cluster mass. 

We therefore conclude that our results are robust to any 
possible aperture effects. 



3 RESULTS 

In the upper panel of Figure[2]we plot the BCG {z < O.I) and 
BGG aperture A'-band luminosities versus their host system 
mass. The galaxies in the lower-mass systems are fainter and, 
morphologically, more likely to be late-type galaxies (Sa to 
Irr) in systems without extended intra-group X-rays. The 
dispersion in the central galaxy luminosities dec reases as the 
host h alo mass increases, as also observed by IZheng et al.l 
|2003). They suggest that this is a result of a broader distri- 
bution of major star formation epochs in lower mass halos. 
This is consistent with our observation that many of the 
galaxies in the low-mass halos are late-types. 

Fitting a relationship to these data we find Lk oc 
M°af^°-'"^ for all BCGs and BCGs. This is steeper than 
has previously been observed. However, given that the ex- 
tended X-ray emission characteristic of the potential well of 
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Figure 1. Upper Panel: Comparing our apparent aperture mag- 
nitudes with the apparent isophotal magnitudes from Lin & Mohr 
(2004). We find an offset of Amu = 0.35 ± 0.03 mag for the 24 
galaxies in common. Lower Panel: Comparing model aperture 
and isophotal luminosities with host system mass. 



a bound system is not detected in the less-massive groups, 
it is more reliable to determine this relationship from the 
early-type BCGs in groups with extended X-ray emission 
and BCGs (eBGG-|-BCG sample). For these we find a rela- 
tionship Lk oc M°nn^^° °^. This is consistent with values 
from the literature (e.g. iLin fc Mohij |2004| : iPopesso et al.l 
l2007i : iHansen et al.ll2007l ). We note that the uncertainties 
in these data have insignificant effect on the fitted rela- 
tions. Analysing the BCGs alone yields a shallower relation: 
Lk oc Mjqq^*'^'^''. This relationship is consistent with no 
increase of BCG luminosity with increasing cluster mass. 
However, due to the scatter in the data it is also consistent 
with the relationship for the eBGG-|-BCG sample and with 
those relationships in the literature. 

Whether the LcMh-relation evolves with redshift is also 
interesting. We examine the 63 BCGs at higher redshift 
(O.I < z < 0.8) in the lower panel of Figure [5] The best- 
fitting relationship to these data is Lk oc 

A'/2oo'*° ''- This 
is entirely consistent with the low-redshift relations. 

In Figure |3] we show how our results compare to those 
measured in the literature, and to the predictions of CMOS 
and VO06. To add anoth er high-redsh if t (z ^ I) point, we 
assume that because the IZheng et al.l l|2007l ) data fit the 
VO06 model, they can be represented by L oc M'^'^^ . There 
is no appreciable offset between the relationship measured 
from HOD analyses or directly from BCG samples. These 
suggest that there is no evolution in this relationship be- 
tween z ~ 0.0 1 and z ^ 1 . Our data are also consistent 
with the predictions of both CMOS and VO06. 
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Figure 2. The relationship of BGG and BCG aperture ii'-band 
luminosities with host system mass. The upper panel shows the 
relationship for the BGGs and the BCGs in low-redshift clus- 
ters (z ^ 0.1). In this panel the points indicate the BCGs (filled 
squares) and the properties of the BGGs (circles denote early- 
type galaxies, triangles denote late-types) and their host groups 
(filled points indicate extended X-ray emission, open points in- 
dicate galaxy halo emission and arrows indicate upper limits 
for the X-ray undetected groups). The solid line indicates the 



0.24±0.08 
200 



relation for the (eBGG+BCG sample), while 



O.liiO.lO 



relation for 



Lk oc M. 

the dot-dashed line indicates the Lk oc M^qq 
the BCGs alone. The lower panel shows the relationship for 
the BCGs in high-redshift clusters (0.1 < z < 0.8). The solid 
line indicates the best-fit Lj( oc Mjqq^^"'^^ relation for these 
high-redshift data and the dashed line indicates the low-redshift 
eBGG+BCG sample for direct comparison. 




Figure 3. Comparing our results with gradients from LcM^- 
relation studies in the literature with redshift. Stars indicate the 
predictions from CMOS and VO06, filled squares indicate previ- 
ous BCG studies, open squares indicate previous HOD studies 
and open circles indicate the results presented here. 



4 DISCUSSION 

BCGs are often thought to be the most conspicuous ex- 
amples of central cluster galaxies undergoing rapid growth 
through merging. These data suggest a picture in which the 
near-infrared luminosity (~stellar mass) of early-type BCGs 
evolves steadily with the growth of their host cluster. In ob- 
serving the LcMh-relation to have a slope ~ 0.3 at higher 
redshift, BCG growth appears to still be in the 'dynamical 
friction time-scale-limited' regime (as per CMOS) at these 
earlier times. Hence, BCG-halo mass growth is controlled by 
this restriction over the cosmic time interval studied here, 
and is not consistent with BC Gs doubling in ma s s from 
z ~ 1 to z ~ as suggested bv lde Lucia fc BlaizotI (|2007l ). 
These observations are more consistent with BCGs accret- 
ing 10 — 20 per cent in stellar mass as their cluster doubles 
in mass. 

This conclusion is consistent with the observation that 
the position of BCGs on scaling relations such as the fun- 
damental plane is a result of having undergone at least 
1 major merger ev ent (e.g. iBovlan-Kolchin et al.l l2006l : 
iBernardi eraLli200'it). These majo r merg er events have been 
directly observed: ^Ya mada et al.l (|2002l ) observed a nearly 
equal-ma ss BCG merger in a ~ 2 x 10^* Mq cluster at 
z ~ 1.26. iRines et al.l (|2007l ) have observed a plume around 
a BCG in a ~ 1 X 10"/i"^Mq cluster at z ~ 0.4. Both 
of these systems lie within the scatter of our data on the 
LcMh-relation. 

Observations of the mass assembly of other massive 
early-type galaxy populations are contradictory. Some claim 
evid ence for no evolution beyond passive for these galaxies 
(e.g. ICimatti et al]|2006l : iBundv et allbood ). whilst others 
clai m stellar mass evolution of a factor of 2 since 2: ~ 1 
(e.g. Ivan Dokkumll2005l : iBell et al.1 [20061 ). In between these 
extremes are studies more consistent with this work, find- 
ing evidence for growth of 50 per cent since 2 ~ 1 
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From a theoretical perspective, we have demonstrated 
that our data are consistent with the predictions of two 
HOD models at z ~ 0.1. However, these models do not 
make predictions for higher redshifts. Examining the evo- 
lution of BCGs in a A -Cold Dark Matter Univ erse, the 
semi-analytical model of lde Lucia fc BlaizotI (|2007l ) predicts 
that ~ 50 per cent of the mass of BCGs is assembled since 
2: ~ 1, with an equivalent growth in the host halo over this 
times cale. This is in c onsist ent with our observations. How- 
ever, lAlmeida et al.l (|2007l ) have analysed the stellar mass 
evolution of large red galaxies between z ~ 0.24 and 2 ~ 0.5 
in the Durham semi-analytic models, and find that, on aver- 
age, large red galaxies assemble ~ 25 per cent of their mass 
since z ^ 1. 

The results presented in this letter are clearly consistent 
with the picture emerging from studies of the evolution of 
massive early-type galaxies and semi-analytic models of the 
evolution of these massive galaxies, that massive galaxies 
have increased in stellar mass by ^ 50 per cent since z ~ 1. 
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5 CONCLUSIONS 

We have examined the BGG and BCG luminosity - host 
halo mass relation at: z < 0.1 and, at higher redshifts than 
previously measured: 0.1 < z < 0.8. We conclude that: 

• BCGs follow a relationship with the mass of their host 
cluster. This relationship does not evolve from z ~ 1 to 
~ 0. 

• Our data are consistent with the predictions of two il- 
lustrative HOD models for the relationship of central galaxy 
luminosity with host halo mass. 

• The BCG - halo-mass relation appears to be controlled 
by the timescale for dynamical friction even at 0.1 < z < 0.8. 
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